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ABSTRACT 

m 
> 

f-^ \ Near-future SZ surveys such as ACT, SPT, APEX, and Planck will soon 

find thousands of galaxy clusters. Multi-frequency arcminute-resolution SZ ob- 
servations can, in principle, determine each cluster's gas temperature (T e ), bulk 
velocity (v pec ), and optical depth (r). However, the frequency bands and detec- 
tor sensitivity employed by upcoming surveys will generally not be sufficient to 
disentangle the degeneracy between these three cluster parameters, even in the 
absence of SZ signal contamination from point sources and imperfect primary 
microwave background subtraction. Assuming contaminants can be removed, we 
find that near-future SZ surveys will be able to constrain well two cluster gas 
parameters that are linear combinations of rT e , tv z , and rT e 2 . Because the SZ 
intensity shift is nearly a linear function of rT e , rv z , and tT £ 2 , a correspondence 
exists between the two effective gas parameters that SZ surveys can constrain 
and simple line-of-sight integrals through the three dimensional cluster. We il- 
lustrate the parameter constraints and correspondence to line-of-sight integrals 
using three dimensional Nbody + hydro cluster simulations and a Markov chain 
Monte Carlo method. We show that adding an independent temperature mea- 
surement to upcoming SZ data breaks the parameter degeneracy and that the 
cluster effective velocity thus constrained is approximately the optical-depth- 
weighted velocity integrated along the cluster line of sight. A temperature prior 
with an error as large as 2 keV still gives bulk velocity errors of 100 km/sec or 
less, even for a more typical cluster with an electron temperature of 3 keV, for 
ACT-like SZ observations in the absence of signal contamination. The Markov 
chain constraints on t> pec and r that we obtain are more encouraging and most 
likely more accurate than those obtained from Fisher matrices. 
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Subject headings: cosmology: theory — large-scale structure of the universe 
- cosmic microwave background - - galaxies: clusters: general - - Sunyaev- 
Zel'dovich effect — parameter extraction 

1. INTRODUCTION 

Observations of the Sunyaev-Zel'dovich (SZ) effect (Sunyaev & Zel'dovich 1972) offer the 
hope of revealing much about the properties of galaxy clusters and the evolution of large-scale 
structure. Several instruments are being built (ACT, SPT, APEX, Planck) that will make use 
of the SZ effect to produce deep galaxy cluster surveys, and upgrades to current experiments, 
such as SuZIE III, will produce deep targeted observations of known galaxy clusters. These 
surveys, in addition to measuring the number density of clusters, can, in principle, reveal each 
galaxy cluster's peculiar velocity, gas temperature, and optical depth, if the SZ information is 
fully exploited. The aim of this paper is to quantify some of the difficulties with determining 
these individual cluster parameters from future SZ measurements and to discuss what cluster 
parameters these future surveys can constrain. 

The SZ effect is a spectral distortion of the cosmic microwave background caused by an 
intervening galaxy cluster. The hot gas in the intracluster medium inverse Compton scatters 
the microwave photons creating this distortion. For reviews of the SZ effect see Rephaeli 
(1995), Sunyaev & Zel'dovich (1980), Birkinshaw (1999), and Carlstrom et al. (2002). The 
dependence of this spectral distortion on the gas temperature (T), radial peculiar velocity 
(v), and optical depth (r), including relativistic corrections, has been computed in several 
papers (Challinor & Lasenby 1998; Sazonov & Sunyaev 1998; Itoh et al. 1998; Nozawa et al. 
1998; Molnar & Birkinshaw 1999; Dolgov et al. 2001). Since the amplitude of this distortion is 
also a function of frequency, SZ measurements at three different observing frequencies would 
ideally be enough to disentangle the three unknowns (T, v, r) contributing to the SZ signal. 
Some first attempts to determine cluster temperatures using the relativistic SZ effect were 
carried out by Pointecouteau et al. (1998) and Hansen et al. (2002), the former simulating 
clusters observed with Planck and the latter providing a temperature estimate of Abell 2163 
with quite large uncertainties. More recently, Hansen (2004a) has developed code to extract 
cluster parameters using the SZ effect which results in rather large errors if given prospective 
data from upcoming surveys. The difficulty, as pointed out in Holder (2004) and Aghanim 
et al. (2003), is that there exist degeneracies among the physically interesting parameters 
(T, v, r) that can best be broken by choosing one frequency relatively low (around 30 GHz) 
and placing the other two around 150 GHz and 300 GHz. In addition, the observations need 
to be of arcminute-resolution to resolve individual clusters. A 30 GHz observing band is not 
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a realistic option for upcoming arcminute-resolution bolometer based instruments (such as 
ACT, SPT, and APEX) since high- sensitivity bolometers sharply lose sensitivity below 90 
GHz and the single-dish diameter required for arcminute-resolution observations at 30 GHz 
is unrealistically large. In addition, large interferometers capable of observing at 30 GHz 
with arcminute-resolution can only view small areas of the sky at a time, making a large 
survey impractical. Thus it is important to determine what information about individual 
galaxy clusters these future SZ surveys will be able to constrain given the reality that they 
will have approximately arcminute-resolution observations at or above 90 GHz. 

In this paper we make a preliminary investigation of information we can potentially 
extract from SZ measurements. In particular, we do not include the effects of a variety 
of real world complications including radio and infrared point sources, galactic dust, non- 
SZ microwave fluctuations or instrumental systematic effects. While these are certainly 
important considerations for any data set, the goal of this paper is to elucidate the maximum 
information about individual clusters we can in principle obtain. Rough estimates of the 
impact of these various contributions are briefly discussed in §8. Detailed analysis including 
these effects is in progress. 

In the next section we briefly summarize the SZ effect. In §3 we investigate how varying 
observing frequencies and detector sensitivity affects parameter degeneracies and parameter 
extraction. In §4 a Markov chain/Fisher matrix method is used to determine which cluster 
parameters can be well constrained by future SZ measurements, and in §5 this method 
is applied to simulated Nbody+gasdynamics clusters and the results presented. In §6 we 
discuss the near-linearity of the SZ intensity shift with respect to rT e , tv z , and rT e 2 , and show 
the resulting close correspondence between the constrained effective parameters from 2D SZ 
images and line-of-sight integrals through the 3D cluster. In §7 we show that an independent 
measurement of T e breaks the parameter degeneracy and use a simple analytical model to 
show that the velocity thus determined is approximately the optical-depth-weighted velocity 
integrated along the cluster line of sight. We then use a Markov chain to calculate errors 
on cluster velocities and optical depths given X-ray temperature priors. In §8 we discuss 
sources of contamination to the SZ signal, and we conclude with a summary of the above. 



2. SZ EFFECT 

When microwave photons pass through the hot gas in the intracluster medium of a 
galaxy cluster, roughly 1% of the photons interact with the free electrons in the gas. These 
photons are inverse Compton scattered and energy is transferred from the hot electrons to 
the cool photons, causing a slight distortion of the microwave background spectrum. This up- 
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scattering of photons causes the intensity of photons with frequencies below about 220 GHz 
to decrease while the intensity of photons with higher frequencies increases. This process is 
called the thermal SZ effect, and it causes an effective temperature shift relative to the mean 
microwave background temperature on the order of one part in 10 4 . If the galaxy cluster 
has some bulk velocity with respect to the microwave background rest frame, then this will 
Doppler shift the scattered microwave photons and cause an additional spectral distortion. 
This further shift in the microwave spectrum is referred to as the kinematic SZ effect and is 
typically an order of magnitude smaller than the thermal SZ effect. The derivation of the 
combined SZ effect can be found in Sunyaev & Zel'dovich (1970) and Sunyaev & Zel'dovich 
(1972), and more recent papers (e.g. Itoh et al. (1998); Nozawa et al. (1998)) have included 
relativistic corrections to these derivations. 

The expression for the SZ effect we use throughout this work is from Nozawa et al. 
(1998) and is given by 
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where I = 2(k B T CMB ) 3 /(hc) 2 , X = hu/k B T CM B, e = k B T e /m e c 2 , v tot is the peculiar 
velocity, v is the peculiar line-of-sight velocity, r = ut J n e dl is the optical depth, and the 
F's and C's are numbers that depend on frequency in a well-defined way. Note that r9 e is 
equivalent to the well known Compton-y parameter which we will show is tightly constrained. 
The most dominant terms in this expression are proportional to rT e , rT e 2 , and tv. It is also 
important to note that the above expression is independent of redshift. This makes the SZ 
effect a powerful probe of the high-redshift universe because the amplitude of the SZ signal 
does not weaken at high redshift (for fixed T, v, r), unlike X-ray and optical signals. Since 
this microwave intensity shift for a given frequency is a non-trivial function of the cluster's 
gas temperature, peculiar velocity, and optical depth, one could hope that choosing at least 
three well placed observing frequencies would allow these three cluster parameters to be 
separated and measured. In practice, the observing frequencies and sensitivities available 
to upcoming SZ surveys result in degeneracies among these parameters. These degeneracies 
are illustrated in the next section. 
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3. PARAMETER EXTRACTION 

3.1. Creating Likelihood Surfaces 

To understand the intrinsic limitations in determining a galaxy cluster's gas tempera- 
ture, peculiar velocity, and optical depth from multi-frequency SZ measurements, we create 
likelihood surfaces for these parameters and compare the 1-a regions for various choices of 
observing frequencies and detector sensitivity. We do this by first assuming some region of 
gas with uniform temperature T e , peculiar line-of-sight velocity v, and total optical depth r. 
We neglect any transverse peculiar velocity since the full SZ effect generates a temperature 
shift on the order of one part in 10 4 , and the transverse velocity component of the SZ effect 
contributes a temperature shift on the order of one part in 10 7 . 

Using the expression for the intensity shift given in eq. (1), we calculate the change 
in intensity one would measure from our fiducial gas region at three different observing 
frequencies. We perform these calculations for several frequency sets. The first observing 
frequency set we choose to be (30, 150, 300 GHz), which is the optimal frequency set found by 
Holder (2004) and roughly that found by Aghanim et al. (2003). We also choose the sets (90, 
150, 300 GHz) and (145, 225, 265 GHz), the latter being the frequencies planned for ACT 
(see Kosowsky (2003)). After obtaining the intensity shift for each of the three frequencies 
in each set, we change variables to the ratios x\ = AI Vl / AI V2 and = AI Vz / AI V2 , where v\, 
and z/ 3 are the three frequencies in a set. These ratios are independent of r (though their 
errors depend on r), so we have two equations and two unknowns (T e and v). We eliminated 
t so that we could more easily construct and view the likelihood surfaces of the remaining 
parameters. Assuming first lfiK detector noise per arcminute beam and then 10 fiK noise 
per beam, we step through the T e — v parameter space and calculate the likelihood of the 
model described by (T e , v), given the underlying fiducial model, in the usual way. Thus, 
we create a 2-dimensional likelihood surface for the parameters T e and v for each observing 
frequency set at both detector noise levels. A 1-a contour is then drawn for each likelihood 
surface by connecting all the points with x 2 — 2.3 (e.g. Press et al. (1997)). 



3.2. Degeneracy Between Cluster Gas Parameters 

The 1-a contours we obtain from these likelihood surfaces verify that the set of frequen- 
cies 30, 150, 300 GHz puts tighter constraints on T e and v than the other sets. With 1/iK 
detector noise per arcminute beam, the set 145, 225, 265 GHz exhibits a clear degeneracy 
between the T e and v parameters, and with 10 fiK detector noise per arcminute beam all 
three frequency sets show significant degeneracy. Figures la-lc show the comparison of 1-a 
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contours for 10 keV (1.2 x 10 8 K) gas with a line-of-sight velocity of 200 km/sec, optical 
depth of 0.012, and 1/iK detector noise. We choose gas parameters corresponding to a hot 
cluster because it is the case which has the highest signal-to-noise. These parameter values 
roughly correspond to the known cluster MS 0451 (Reese et al. 2002). The gas temperature 
is constrained to within 0.5 keV and its velocity to within 25 km/sec at the 1-a level for 
30, 150, and 300 GHz observing frequencies (figure la). Shifting the lowest frequency to 90 
GHz increases these uncertainties by a factor of two (figure lb). Using observing frequencies 
at 145, 225, and 265 GHz results in temperature uncertainties of 6 keV and velocity un- 
certainties of 220 km/sec (figure lc). Figures 2a-2c show the 1-a likelihood contours for 10 
keV gas with a line-of-sight velocity equal to -200 km/sec, optical depth of 0.012, and 1/iK 
detector noise. Comparing these to the previous group of figures indicates that a negative 
line-of-sight velocity somewhat increases the uncertainties obtained for all frequency sets, 
as has been pointed out in Aghanim et al. (2003). Adding a 30 GHz or 90 GHz observing 
frequency to the 145, 225, 265 GHz frequency set with lfiK noise greatly reduces the 1-a 
regions and makes the temperature and velocity constraints as tight as for the 30, 150, 300 
GHz and 90, 150, 300 GHz frequency sets respectively. These results confirm previous results 
by Holder (2004) and Aghanim et al. (2003). It is clear that measurements near the null 
of the SZ effect are not particularly useful for SZ studies from a signal-to-noise perspective, 
but such measurements will be useful for microwave observations that aim to minimize SZ 
contamination and as a useful diagnostic of point source contamination. 

With 10/iK detector noise per arcminute beam, figures 3a-3c show the 1-a contours for 
the same 10 keV gas region with 200 km/sec line-of-sight velocity and an optical depth of 
0.012. The gas temperature is constrained to within 5 keV and the velocity to within 200 
km/sec at the 1-a level for the 30, 150, and 300 GHz set (figure 3a). Increasing the lower 
frequency to 90 GHz constrains the temperature to within 8 keV and the velocity to within 
250 km/sec (figure 3b). The 145, 225, and 265 GHz frequency set gives 1-a errors of 9 keV 
and 350 km/sec for temperature and velocity respectively with 10 fiK detector noise (figure 
3c). Thus to constrain the cluster parameters well using SZ observations alone (in the absence 
of SZ signal contamination from point sources or residual primary microwave background), 
one needs both a low frequency band (90 GHz or lower) and a detector sensitivity not much 
higher than 1/iK per arcminute beam. 

We also compare the 1-a likelihood contours for gas with different Compton-y parame- 
ters keeping the observing frequencies fixed. In figure 4a we see the 1-a contour for a 7 keV 
gas region with an optical depth of 0.009, 200 km/sec line-of-sight velocity, 1/iK detector 
noise, and observing frequencies at 30, 150, and 300 GHz. Figure 4b shows the 1-a contour 
for a 3 keV gas region with an optical depth of 0.004 and the same line-of-sight velocity, 
detector noise, and observing frequencies. The 1-a contour for the gas region with lower 
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Compton-y parameter shows a larger parameter degeneracy for the same set of observing 
frequencies and detector noise. This is because a decrease in Compton-y parameter lowers 
the signal-to-noise. The signal-to-noise is reduced by two effects as one moves to lower mass 
clusters: a lower signal in all bands and reduced relative importance of relativistic effects. 
We have verified that the former is a more important contributor to parameter degenera- 
cies than the latter. Table 1 lists the above 1-a errors on T e and v for differing observing 
frequencies, detector noise, and Compton-y parameters for convenient reference. 

4. CONSTRAINED PARAMETERS FROM FUTURE SZ SURVEYS 

Despite the inability of upcoming experiments like ACT, SPT, and APEX to determine 
all of the cluster gas parameters, they will provide tight constraints on certain combinations 
of parameters. The next two sections explicitly demonstrate the parameter combinations 
which will be determined with good precision by these kinds of experiments. 

4.1. Markov Chain Analysis 

We create a Markov chain Monte Carlo (MCMC) using the parameters T e , v, and 
t to find realistic error regions for all three parameters. This MCMC is made using the 
Metropolis-Hastings algorithm which randomly steps through a parameter space and accepts 
all points whose likelihood is greater than the previous point. If the likelihood is less than 
the previous point, the current point is accepted with a probability given by the ratio of the 
two likelihoods. A comprehensive review of MCMCs can be found in Gilks et al. (1996); 
they were introduced into cosmology by Christensen & Meyer (2001), Christensen et al. 
(2001), Lewis & Bridle (2002), and Kosowsky et al. (2002). The region in the 3-dimensional 
parameter space containing 68% of all the points accepted in the chain we define as the 1-a 
region. This 1-a region is projected onto 2 dimensions in the figures below. Figures 5a-5c 
show the 1-a regions generated by a MCMC for 10 keV gas with 200 km/sec line-of-sight 
velocity and an optical depth of 0.012. The frequency set 145, 225, and 265 GHz was used 
with lfiK detector noise per beam to simulate ACT observations. In the Markov chain, the 
parameter space was restricted to T e e (0,2 x 10 8 K (17 keV)), r e(0, 0.02), and v €(-1500 
km/sec, 1500 km/sec), and 5 million steps were used with about 50% acceptance rate in 
the chain. In addition, each step was taken in all three parameter directions simultaneously 
with different step sizes in each parameter direction. From the 1-a regions in figures 5a-5c, 
we can see a clear degeneracy among all three parameters. 
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Figures 6a-6c show the same 1-a region as figures 5a-5c except using the parameter 
directions rT e , rT e 2 , and tv. We have verified that the MCMC in either variable set gives 
the same results, although it is significantly more efficient to use the second set of variables. 
We choose this set of parameters since these are the dominant terms in the intensity shift 
expression in eq. (1). These figures demonstrate that the cluster parameters lie in a nearly 
1-dimensional subspace of the 3-dimensional parameter space given by rT e , rT e 2 , and tv. We 
can find two axes within this parameter space (both orthogonal to the degeneracy direction) 
along which the cluster parameters are tightly constrained and one axis (parallel to the 
degeneracy direction) along which the cluster parameters are largely unconstrained. 



4.2. Fisher Matrix Determination of Constrained Parameters 

The three orthogonal directions in this parameter space that allow us to tightly constrain 
the gas parameters in two directions with one direction unconstrained correspond to the 
principal axes of the 1-a error ellipsoid which is projected in 2 dimensions in figures 6a-6c. 
To find these principal axes we calculate a Fisher matrix at the most likely point found by 
the MCMC. The Fisher matrix describes the curvature of the likelihood surface at a given 
point in parameter space. It is expressed by the formula 



Since we have C oc e x2 / 2 , the Fisher matrix becomes F a p = \ , d x , . Using 
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where p = (rT e /(10 6 K), rv/(l km/sec), rT 2 /(10 14 K 2 )), AI°. is the observed AJ„., and e„. 
is the error on A7°., and averaging gives 
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(See Dodelson (2003) for a good explanation of the Fisher matrix and its applications.) The 
eigenvectors of the Fisher matrix evaluated at the minimum point of x 2 correspond to the 
principal axes of the error ellipsoid. 

This technique is essentially a version of principal component analysis (PCA). In PCA 
it is known that the principal components one finds depend on how the variables are scaled. 
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There is no single correct scaling to choose, but the one that is widely preferred in physical 
uses is scaling all the variables to order unity, which has the advantage that all the variables 
are weighted the same. (Jackson (1991) provides a good overview of PCA.) We scale rT e , 
tv, and tT 2 by 10 6 K, 1 km/sec, and 10 14 K 2 respectively, which are characteristic cluster 
values for these variables. 

The principal axes are linear combinations of the previous parameter directions such 

that 

a = C(p,)p, (5) 

where the rows of C(pJ are the eigenvectors of F, a = (a, b, c) is a point in the new parameter 
space, and p are the old parameter directions described above. The vector p^ describes p 
evaluated at the maximum likelihood point. Note C(pJ (and thus a) will differ with p^, 
but the variation with p^ is fairly weak for realistic parameter regions. Roughly, the Fisher 
matrix eigenvectors are e x = (1,0,0), e 2 = (0,0.4,0.9), and e 3 = (0,0.9,0.4) normalized 
to unity. Therefore the a parameter is dominated by rT e , and the b and c parameters 
are primarily linear combinations of tv and rT e 2 . Thus, SZ measurements provide precise 
measures of the y parameter and one linear combination of the kinematic SZ effect and the 
relativistic corrections. 

Figures 7a-7c show that the gas parameters for this fiducial model are constrained to 
within 1% in the a direction, 3% in the b direction, and 70% in the c direction. For a 4 
keV gas region with 200 km/sec line-of-sight velocity, an optical depth of 0.005, and detector 
noise of IfiK, the gas parameters are constrained to within 4% in the a direction, 16% in 
the b direction, and 200% in the c direction. The gas parameters are constrained to within 
8%, 22%, and 160% in the a, b, and c directions respectively for a 10 keV gas region with 
200 km/sec line-of-sight velocity, an optical depth of 0.012, and 10fiK detector noise. 

The a and b parameters, which are linear combinations of rT e , tT 2 , and tv, are therefore 
well constrained by an ACT-like SZ survey with observing frequencies near 145, 225, and 265 
GHz and IfiK detector noise. This technique can be applied to determine the constrained gas 
parameters from any multi-frequency SZ observations that are without arcminute-resolution 
observations at frequencies below 100 GHz or that have low frequency information but with 
a noise component considerably larger than 1/iK. By combining this gas information from 
SZ surveys with a data set that can constrain just one of the parameters T e , r, v (such as an 
X-ray survey of clusters that constrains T e for each cluster), all three cluster gas parameters 
can be well determined. 
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5. RESULTS USING SIMULATED CLUSTERS 

We now apply this technique to simulated SZ observations which we generate using 
simulated galaxy clusters. In §5.1 we describe the two simulated clusters we use, one about 
9 keV and the other about 3 keV, and in §5.2 we discuss how SZ simulations are created 
from these. The results of applying the technique in §4 to simulated ACT-like SZ maps of 
both clusters and to simulated Planck-like SZ maps of the 9 keV cluster are contained in 
§5.3. 

5.1. Cluster Simulations 

The clusters we use are two high-resolution 3D cluster simulations that were made using 
the Adaptive Refinement Tree (ART) Nbody+gasdynamics code (Kravtsov 1999; Kravtsov 
et al. 2002). These clusters were simulated using a ACDM model with fl m = 0.3, Vl b = 0.043, 
h = 0.7, and <r 8 = 0.9. Each cluster has a redshift of z=0.43 and is contained in a cube 
of side length 2 Mpc (~ 6 arcminutes). Each grid element within the larger cube has a 
side length of 0.0078 Mpc (~ 0.02 arcminute). The simulations track the density of dark 
matter particles, the density of gas particles, the gas temperature, and the 3-dimensional 
gas velocity for each grid element. One cluster is similar in size to the Coma cluster and has 
a mass of ~ 1O 15 M , an optical depth of ~ 0.01, and an average gas temperature of 9 keV. 
The other cluster is similar in size to the Virgo cluster and has a mass of ~ 2 x 1O 14 M , 
an optical depth of ~ 0.005, and an average gas temperature of 3 keV. Both clusters have 
characteristic bulk velocities of several hundred km/sec. The morphology of the Virgo-size 
cluster indicates that it consists of a recent merger of two smaller clusters. These cluster 
simulations do not include the effects of gas cooling, stellar feedback, magnetic fields, and 
thermal conduction. For a more detailed cluster description see Nagai & Kravtsov (2003) 
and Nagai et al. (2003). 

5.2. SZ Map Generation 

To create a simulation of an SZ observation, we choose one of the cluster simulations and 
a set of observing frequencies. The SZ intensity shift that microwave photons would experi- 
ence passing through each grid element is calculated using eq. (1) and the gas temperature, 
gas density, and gas velocity of each element. Every Al is then integrated over a frequency 
band centered around each observing frequency in the set. We use a 3 GHz frequency band- 
width, as opposed to a more realistic bandwidth of ~ 25 GHz, for numerical convenience. 
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However, the bandwidth size has a negligible effect on cluster constraints, which we verified 
by redoing some of our results with a 25 GHz bandwidth. Thus, for each frequency band we 
end up with an SZ cube of AI values. This cube is then projected along the line of sight into 
a two-dimensional SZ distortion of the sky. We do not include the primary microwave back- 
ground in our simulations, assuming it is perfectly subtracted, since it varies on scales large 
compared to the cluster. The main effect of residual primary microwave contamination will 
be as a source of noise for extracting estimates of the peculiar velocity from the constrained 
parameters. Figures 8a-8c and 9a-9c show the 2D SZ images of the 9 keV and 3 keV clusters 
after this projection process for the frequency bands centered on 145, 225, and 265 GHz. 
After creating a 2D SZ image for each observing frequency band, we smooth each image 
by convolving it with a Gaussian beam, and increase the pixel size of our images to sizes 
realistic for upcoming SZ surveys by averaging together smaller pixels. Finally, Gaussian 
random noise of standard deviation equal to our chosen detector sensitivity is added to each 
2D pixel. 

We make SZ simulations using ACT instrument specifications for both the 3 keV and 
9 keV clusters and using Planck specifications for the 9 keV cluster only. The ACT-like 
SZ images of the 9 keV cluster (assuming perfect microwave background and point source 
removal) are shown in figures 10a- 10c. For these images we use the frequency bands centered 
on 145, 225, and 265 GHz. We assume a beam size of 1 arcminute and choose a pixel size of 
0.3' x 0.3'. Each 0.3 arcminute pixel is given 3fiK of Gaussian random detector noise. Our 
results are not qualitatively different for moderately different noise levels between channels. 
Figures lla-llc show SZ simulations of the 3 keV cluster with ACT specifications as above. 
For comparison, we also made maps appropriate to the Planck experiment, with a beam 
size of 4 arcminutes and detector noise of 16fiK per 2' x 2' pixel, at the same frequencies as 
the ACT maps (which are similar to the actual Planck bands centered at 143, 217, and 353 
GHz). 



5.3. Parameter Constraints from Simulated 3D Clusters 

We now apply the Markov chain/Fisher matrix technique described in §4 to the simu- 
lated ACT-like and Planck-like SZ images of the simulated clusters. The only change is that 
in eq. (3) we assume 3/j.K noise for ACT-like images and !6fiK noise for Planck-like images. 
The cluster parameters we constrain using this method are really effective parameters that 
correspond to integrals of the cluster parameters along a line of sight. The SZ intensity (eq. 
[1]) is not linear, so there is no guarantee that the sum of SZ signals of varying temperature 
and velocity can be fit to a single temperature and velocity. We fit the resulting SZ intensity 
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as a function of frequency to a model with a single temperature and velocity and call the 
constrained parameters a e s; we discuss what integrals these effective parameters correspond 
to within the three dimensional cluster in §6. 

Figures 12a-12c show the projected 1-a contours for the a cS parameters for the central 
pixel of ACT-like SZ images of the simulated 9 keV cluster. We do this analysis on a pixel- 
by-pixel basis to potentially obtain the most information about cluster substructure. (See 
Nagai et al. (2003) for some discussion of substructure.) In a low signal-to- noise experiment 
or as a means to average out substructure, it could be advantageous to add together pixels. 
The 1-a errors on a e g, b e s, and c e fr for this central pixel are 0.06, 0.5, and 5 respectively. The 
range of parameters T G (0, 2 x 10 s K), r G (0,0.02), and v G (-1500 km/sec, 1500 km/sec) 
correspond to ranges of about a G (0,4), b G (—5,20), and c G (—30,24). Clusters with 
larger SZ signals (and thus larger T e 's and r's) tend to have larger a, b, and c parameters, 
and this is borne out in the comparison of our results for the 9 keV and 3 keV clusters. 
Since the a, b, and c parameters have been scaled to roughly order unity for characteristic 
cluster values, the absolute errors are a meaningful reflection of how well the gas properties 
are constrained. 

Similar results are obtained from the ACT-like SZ simulation of the 3 keV cluster and 
the Planck-like SZ simulation of the 9 keV cluster. Figures 13a-13c show the projected 1-a 
contours for the central pixel of the ACT-like SZ images of the 3 keV cluster. These figures 
demonstrate a aeS ~ 0.02, o& cff ~ 0.2, and a CcS ~ 1.5. Thus a e s, b e g, and c c g are constrained 
to a small region of the available parameter space, with the first two components especially 
well-constrained. Figures 14a-14c show the projected 1-a contours for the central pixel of the 
Planck-like SZ images of the 9 keV cluster. These figures show a acH ~ 0.06, a bcH ~ 0.7, and 
a CeS ~ 3, indicating again that the measurements are providing strong constraints within 
the available parameter space. 

These results demonstrate that the a e fr and b e g cluster parameters are well constrained 
and c cS is moderately well constrained by SZ observations typical of ACT and Planck, 
assuming perfect primary microwave background and point source removal. 

6. CORRESPONDENCE BETWEEN CONSTRAINED EFFECTIVE 
PARAMETERS AND LINE-OF-SIGHT INTEGRALS 

The cluster gas parameters we have constrained in the previous sections using a pro- 
jected two dimensional SZ image correspond to integrals along the line of sight of the three 
dimensional cluster. Previous studies of cluster projection effects have been done indicating 
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that the temperature in principle obtainable from SZ measurements is really a Compton- 
averaged quantity (i.e. each line-of-sight integral of T e is weighted by the Compton pa- 
rameter) (Hansen 2004b; Knox et al. 2004). Here we show that the line-of-sight integrals 
corresponding to the a eS parameters are even more straightforward. 

The reason for the simple correspondence is that the SZ intensity shift given by eq. (1) 
is nearly linear with respect to tT, tv, and tT 2 . These three terms are the most dominant 
terms in the expression and represent most of the change in intensity for temperatures of 
several keV and velocities of several hundred km/s. If AI U were exactly a function only of 
tT, tv, and rT 2 in eq. (1), then AI V would be exactly a linear function of a, b, and c. In 
that case, the measured a e s would be equal to Cpj, where C is given in eq. (5), and the 
sum is over the gas properties p i of each element i along the line of sight. 

The full SZ intensity shift expression in eq. (1) includes terms non-linear in tT, tv, and 
tT 2 . However, the addition of these non-linear terms in the calculation of AI V integrated 
along a line of sight only introduces a slight bias between a cf[ and J Cdp t . In figures 12-14, 
within the 1-er contours of a cS , we indicate the best fit a eff found by the MCMC, for a 
given noise realization, by a star (-fc). The best fit a cff obtained if the clusters are observed 
with an ideal instrument without detector noise are indicated by dots (•). Diamond shapes 
(♦) mark the values of the line-of-sight integrals given by J Gdp i and calculated using the 
3D cluster simulations. In figure 12a, the difference between a e s from an ideal, no noise 
instrument and from the line-of-sight integral is Aa = 0.001. The difference between b c g 
from an ideal instrument and from the line-of-sight integral is Ab = —0.03. Clearly, in 
the absence of detector noise, the correspondence between a c s and b c s and the line-of-sight 
integrals is very close. Moreover, the difference between a e s and 6 e fr and the line-of-sight 
integrals given realistic detector noise is still well within the 1-er errors on a c g- and b cS 
given by the MCMC. Similar results can be seen in figures 13a and 14a. These simulations 
demonstrate an agreement to within 1-a between the a c s and b e s parameters constrained by 
SZ measurements and line-of-sight integrals given by f Cdp { . 

Figure 15 shows this correspondence more explicitly. Plotted on the y-axis are best-fit 
a cS obtained via a MCMC using intensities from simulated noise-free SZ images of the 9 
keV and 3 keV clusters assuming ACT-like observations. On the x-axis are plotted J Gdp { 
for the corresponding lines of sight. Four different lines of sight through both the 9 keV 
and 3 keV simulated clusters are plotted. These lines of sight are 0', 1', 1.5', and 2' from 
the central pixel of the simulated SZ images. Filled shapes correspond to the 9 keV cluster 
and unfilled shapes to the 3 keV cluster. For the a c s and b cS parameters, which can be well 
constrained, their equivalence to J Qdp i is extremely close. This confirms the near linearity 
of the SZ intensity shift with respect to the a, b, and c parameters. The c e ff parameters 
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demonstrate some scatter around the y = x line, and this is because the degeneracy in the 
c direction prevents a MCMC from settling on the correct c cff value. 

If we could tightly constrain a e s, b e s, and c c g- via SZ measurements, we could solve for 
the quantities tT t , tv t , and r(T 2 ) r , where the subscript r corresponds to optical-depth- 
weighted integrals (e.g. T T = jTdr/ J dr). From these one can find T y , rT T /T y , and 
v T T y /T T , following the algebra in Knox et al. (2004), where the subscript y corresponds 
to a pressure- weighted integral. Therefore SZ measurements would constrain the pressure- 
weighted temperature (arising from the relativistic corrections), the optical-depth- weighted 
velocity times a correction factor and the optical depth times a similar correction factor 
which is the ratio of different weighted temperatures. Since degeneracies will allow SZ 
measurements to constrain only a c s and b e g, information from an external source will be 
needed to constrain the above physically interesting quantities. 



7. BREAKING PARAMETER DEGENERACY WITH AN X-RAY 

MEASUREMENT OF T e 

7.1. Measured Effective Velocity is approximately J vdrj J dr 

Assuming contamination sources can be dealt with effectively, future SZ observations 
should be able to constrain two quantities given by 

a e ff « r(ciT T + c 2 (T 2 ) T ) + c 3 rv T and (6) 
b cS « t{c±T t + c 5 {T 2 ) t ) +c 6 tv t , (7) 

where the c's are elements of the matrix C defined in eq. (5). It is conceivable that tem- 
perature measurements from an X-ray survey of clusters could allow the determination of r 
and v T . 

Formally X-ray observations would need to provide a constraint on CiT T + Cj(T 2 ) T , where 
the c's are known constants. However, T T and (T 2 ) T are not obviously given by X-ray 
observations. An X-ray derived temperature is also not equivalent to T T . The two may differ 
by as much as 1 keV (Mathiesen & Evrard 2001). If X-ray observations gave T x = T T and it 
was true that (T T ) 2 = (T 2 ) T , then the effective velocity we would get from a MCMC, after 
adding the T x prior, would be equal to v T . For our 9 keV simulated cluster, (T T ) 2 and (T 2 ) T 
agree to within 5% on average, and for our 3 keV simulated cluster, the two agree to within 
12% on average. 

To get an estimate of the biases incurred by not having the correct weighted tempera- 
tures, we assume (T T ) 2 = (T 2 ) T and add a temperature prior of T T to our MCMC. Adding a 
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T T prior, in the manner we discuss further in §7.2, we find for the central pixel of the 9 keV 
cluster, from simulated, noise-free, ACT-like SZ images, an effective velocity of 230 km/sec 
from the MCMC and an optical-depth- weighted line-of-sight velocity of 218 km/sec from the 
three dimensional cluster simulation. For the central pixel of the 3 keV cluster SZ image, 
we find an effective velocity of -10 km/sec from the MCMC and an optical-depth-weighted 
line-of-sight velocity of -3 km/sec from the three dimensional cluster simulation. The bias 
between the velocity from the MCMC and v T is most likely due to the breakdown of the 
(T T ) 2 = (T 2 ) T assumption. 

To quantify the bias incurred from T x differing from T T , we add ail keV offset to our 
T T prior. We find for the central pixel of the 9 keV cluster SZ image, an effective velocity of 
228 km/sec from the MCMC for both + and - 1 keV offsets. We find for the central pixel 
of the 3 keV cluster SZ image, an effective velocity of -45 km/sec from the MCMC for a + 
IkeV offset and an effective velocity of 6 km/sec for a - 1 keV offset. This would suggest a 
total bias between the measured effective velocity and v T of about 15 km/sec for the 9 keV 
cluster and between 10 and 40 km/sec for the 3 keV cluster. 

7.2. MCMC Errors on v c s and r Given an X-ray T cff Prior 

To determine realistic errors on t> cff and r from adding a measurement of T eff to our 
simulated SZ data, we again use a MCMC. We weight each point in our MCMC by the factor 
e -((T-r eff )/Ar eff ) /2 w ] iere we ca i cu l a te T cS = J T e dr / J dr from the 3D cluster simulation and 
AT cS is the assigned measurement error on T cS . 

For a 1 keV error on T eff , our ACT simulation of the 9 keV cluster gives a v = 20 
km/sec and oy = 0.002 for the central pixel. Our ACT simulation of the 3 keV cluster 
gives a v = 60 km/sec and o T = 0.004 for the central pixel. Assuming a 2 keV error on T e s, 
the ACT simulation of the 9 keV cluster gives a v = 40 km/sec and o T = 0.004, and the 
ACT simulation of the 3 keV cluster gives a v = 100 km/sec and o> = 0.005. Our Planck 
simulation of the 9 keV cluster gives a v = 500 km/sec and o T = 0.0024 for a 1 keV error on 
T e ff and a v = 560 km/sec and o T = 0.0062 for a 2 keV error, for the central pixel. Table 2 
lists these 1-a errors for convenient reference. 

The errors obtained on a v and o T from our MCMC are smaller than those obtained 
using a Fisher matrix. However, the MCMC errors are more accurate than those from a 
Fisher matrix since a Fisher matrix approximates the likelihood surface by an ellipsoidal 
Gaussian, which can result in overestimated errors for likelihood surfaces with strong spatial 
curvature such as these. These results show that, in the absence of contamination from 
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imperfect point source and primary microwave background removal, adding X-ray tempera- 
ture measurements to the data from upcoming ACT-like SZ surveys can determine cluster 
peculiar velocities to within 100 km/sec or less. Large-scale velocity fields obtained from 
galaxy clusters out to high redshift could provide an interesting probe of dark matter and 
dark energy. 

8. SOURCES OF CONTAMINATION TO THE SZ SIGNAL 

Major sources of possible contamination have been neglected in this exercise. In par- 
ticular, it is possible that primary and secondary microwave fluctuations and point sources 
(both radio and infrared) could be problematic. 

Primary fluctuations are in some ways both the largest and the smallest concern. The 
fluctuation amplitudes are on the order of 100 fiK and have the exact same spectral behavior 
as the kinematic SZ effect, providing a noise source that is an order of magnitude larger than 
the signal. However, these fluctuations will be highly coherent over the extent of the cluster, 
so the pixel-by-pixel component separation will naturally measure this extended emission. 
At that point, a simple spatial filter can be applied to remove the primary microwave fluctua- 
tions. This spatial filter will have the effect of removing roughly half of the cluster kinematic 
SZ signal (Holder 2004), thereby reducing the signal-to- noise by approximately this same 
factor. 

Secondary fluctuations will be dominated by kinematic SZ from the quasi-linear regime 
(the Vishniac effect; Vishniac (1987)) and the thermal SZ background. The kinematic SZ 
fluctuations are expected to be below the pixel noise, and therefore subdominant, the thermal 
SZ background will simply add noise to the component that is separated as thermal SZ. The 
rms is expected to be roughly 10 ^K, which will most likely serve as the dominant source of 
noise for this component. However, this is much smaller than the expected thermal SZ signal 
from each cluster and should not impact the component separation process at a noticeable 
level. 

Radio point sources that are uncorrelated with galaxy clusters are not a concern (Knox 
et al. 2004), but radio point sources within galaxy clusters can "fill in" the SZ decrement and 
severely impact cluster SZ measurements. This is a long-standing concern for low frequency 
SZ measurements (Moffett & Birkinshaw 1989), and could be a concern even at frequencies 
as high as 150 GHz. The spectra of radio sources up to such high frequencies are not well 
known, but very rough estimates can be made of the most likely contamination. Radio 
surveys at 1.4 GHz have been done of nearby Abell clusters (Ledlow & Owen 1996) and 
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distant X-ray selected clusters (Stocke 1999) that find that a typical galaxy cluster has of 
order one radio source at 1.4 GHz that would be of order a few mJy at a cosmological 
distance. Detailed studies of spectra of bright radio sources indicate (Herbig & Readhead 
1992) that typical radio sources have spectra that are falling and steepening with frequency. 
In particular, in their sample they found only a handful of sources that were as bright at 
40 GHz as at 1.4 GHz. Most sources with rising spectra at 1.4 GHz eventually turned over 
and had lower fluxes at 40 GHz than at 1.4 GHz, indicating that studies based on spectral 
indices at low frequency will not provide accurate estimates of behavior at high frequencies. 
Assuming that cluster sources have the same spectral behavior and that the steepening at 
high frequencies continues, this would indicate that only a few percent of clusters will have 
a radio source contributing of order mJy flux at 150 GHz. In this handful of clusters, radio 
sources will be a concern, as this flux would be an order of magnitude larger than the pixel 
noise. However, in the majority of clusters the contamination due to cluster radio galaxies 
would be comparable to or smaller than the pixel noise. 

Infrared point sources, consisting largely of dusty star forming galaxies, are likely to be 
a major source of contamination (Knox et al. 2004; White & Majumdar 2004). A detailed 
treatment is beyond the scope of this paper, but the results of Knox et al. (2004) suggest 
that hot clusters (with temperatures above about 6 keV) will have a large enough SZ signal 
that the infrared point sources can be estimated simultaneously with the thermal SZ and 
kinematic SZ signals, assuming an independent measure of the gas temperature. In this work 
we have found that the three measurements by an ACT-like experiment provide only two 
effective constraints, suggesting that there is redundancy in the measurements and that an 
additional degree of freedom (infrared point sources) could be allowed without significantly 
degrading the constraints. The surest solution is to use ALMA to measure the relevant point 
source fluxes. As a point of comparison, ALMA could image 100 square degrees (comparable 
to the ACT survey area) to a point source sensitivity of 0.1 mJy at 140 GHz in less than 
1 month. If one were to instead focus on the inner 2' of the largest 100 galaxy clusters, 
this would take several hours of observing time. Note that these same observations could 
be used to estimate the SZ effects, but the small primary beam of ALMA (due to the large 
telescopes) would require careful mosaicing of the cluster to avoid resolving out much of the 
cluster flux. 



9. DISCUSSION AND CONCLUSIONS 



Instruments such as ACT, SPT, APEX, and Planck will find thousands of galaxy clus- 
ters in the near future via SZ observations. In addition to determining the number density of 
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clusters, which can put limits on cosmological parameters, these surveys will reveal informa- 
tion about the gas properties of individual clusters. Ideally, SZ observations would be made 
in at least three frequency bands with one frequency around 300 GHz, one around 150 GHz, 
and another either near 90 GHz or better yet near 30 GHz. Arcminute-resolution observa- 
tions at those frequencies with ~ l/j,K detector noise would tightly constrain the cluster gas 
temperature, line-of-sight velocity, and optical depth in the absence of excessive point source 
and primary microwave background contamination from imperfect subtraction. Without 
this set of SZ observations, parameter degeneracies prevent disentanglement of these three 
cluster parameters. 

Current limitations in technology and instrument availability will make it impractical 
to obtain 30 GHz, arcminute-resolution, 1/iK sensitivity, SZ observations of the majority 
of the clusters that will be found. SZ surveys that will have 90 GHz channels will still 
have parameter degeneracies resulting from detector noise ~ 10/iK. However, we find that 
upcoming SZ surveys will be able to tightly constrain two cluster gas parameters which are 
linear combinations of rT e , tv, and rT e 2 . The constrained parameters are roughly rT e and a 
single linear combination of the other two terms. We demonstrated that this is the case for 
both individual isothermal gas regions and for 3D simulated Nbody + hydro clusters. 

The SZ intensity shift that microwave photons experience passing through a cluster 
is nearly a linear function of rT e , tv, and rT e 2 , these being the most dominant terms in 
the intensity shift expression. This near-linearity results in a close correspondence between 
the two effective parameters SZ surveys will constrain and simple line-of-sight integrals of 
these parameters through the three dimensional cluster. We illustrated this correspondence 
with our three dimensional cluster simulations. This will greatly simplify data analysis of 
multi-frequency SZ data: it will not be necessary (or useful) to model the intensity as a 
superposition of elements along the line of sight but instead the SZ effect can be modeled as 
a single gas element with a single effective temperature and velocity. 

We have shown that a temperature constraint added to SZ data breaks the parameter 
degeneracy between r, T eff , and t> cff . Using the above linearity, we show that the effective 
velocity constrained by combining SZ with an independent temperature measure is approx- 
imately the optical-depth- weighted velocity integrated along the cluster line of sight. Since 
X-ray derived temperatures do not give us precisely the weighted temperature measurements 
that are required to determine J vdT / J drr exactly, we find the measured effective velocity 
will be biased away from f vdr/ f dr by about 15 to 40 km/sec, with a smaller bias for 
hotter, relaxed clusters. 

Errors on r and v e s are calculated via a Markov chain Monte Carlo method assuming 
a temperature prior in addition to SZ data. We find for ACT-like SZ simulations of our 
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9 keV cluster, a v = 20 km/sec and a T = 0.002 for a 1 keV error on T cS , and a v = 40 
km/sec and oy = 0.004 for a 2 keV error on T e g. For our 3 keV simulated cluster, oy = 
60 km/sec and ay = 0.004, and a v = 100 km/sec and oy = 0.005 for 1 keV and 2 keV 
errors on T e g respectively. The Markov chain errors we find on v e g and r are smaller than 
those obtained via a Fisher matrix. A Fisher matrix overestimates the errors because the 
likelihood surface is strongly curved in this parameter representation, strongly violating the 
implicit assumption of ellipsoidal symmetry over the parameter region of interest. Note that 
the errors on velocities will be increased when residual primary microwave contamination 
is included, and that bulk flows within the clusters provide comparable noise in matching 
observed peculiar velocities to the true bulk velocity of the cluster (Nagai et al. 2003; Holder 
2004; Diaferio et al. 2005). 

If an independent cluster temperature estimate from X-ray spectroscopy is unavailable, 
temperature estimates can also be obtained from either the cluster velocity dispersion (Lubin 
& Bahcall 1993) or the cluster integrated SZ flux (Benson et al. 2004) and scaling relations. 
Since an accuracy of only 2 keV is required on an additional temperature measurement to 
obtain very interesting velocity estimates, the use of these scaling relations could prove to 
be a very beneficial tool. 

As discussed in §8, contamination from primary microwave background fluctuations 
and point sources add another source of noise that must be factored into these parameter 
constraints. Radio point sources due to emission from galaxy cluster members themselves 
and infrared point sources will both be non-negligible sources of SZ signal contamination. 
Studies of the effect point source contamination will have on cluster parameter extraction 
have been carried out by Knox et al. (2004) and Aghanim et al. (2004). Both studies have 
found that the contamination could potentially be serious; however the latter study considers 
the effect of point source contamination if no attempt is made to filter point sources out of 
the observations or model them into the parameter extraction routines. Moreover, even in 
the worst case point source contamination scenario, observations with an instrument such as 
ALMA will allow straightforward point source subtraction from SZ images. Clearly either 
filtering techniques or additional ALMA type observations will be needed to minimize both 
the point source and primary microwave background contamination of SZ signals. 

Near-future SZ surveys will open the door to a wealth of information about galaxy 
clusters. Determining the number density of galaxy clusters as a function of redshift is po- 
tentially a strong probe of dark energy's equation of state and variability over time (Haiman 
et al. 2001; Holder et al. 2001; Hu 2003; Majumdar & Mohr 2003, 2004). However, galaxy 
clusters offer more information that can also yield cosmological information. The kinematic 
SZ signature of galaxy clusters can reveal large-scale velocity fields out to high redshift 
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that can provide an alternative probe of large-scale dark matter and dark energy (Peel & 
Knox 2003). Cluster optical depth information yields cluster gas mass estimates, and optical 
depths are crucial to any of the tests that have been proposed using the (extremely difficult 
to measure) polarization of scattered microwave photons at the position of galaxy clusters 
(Kamionkowski & Loeb 1997). The gas parameters T e , r, and v of individual galaxy clusters 
are of direct interest for cluster astrophysics. Arcminute- resolution SZ observations can be- 
gin to probe cluster substructure and offer more information about cluster gas profiles and 
internal gas dynamics. In summary, SZ observations are entering new territory, where large 
scale surveys will offer new understandings of galaxy clusters and cosmology. 

The authors would like to thank Daisuke Nagai for generously providing us with his 
cluster simulations. NS would also like to thank Rouven Essig for useful discussions. This 
work was supported by the NSF award AST-0408698 for the ACT project. 
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Table 1: The 1-a errors on T e and v for different observing frequencies, detector noise, and 
gas parameters. 
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Table 2: The 1-a errors on r and v for ACT-like and Planck-like experiments using 9 keV 
and 3 keV simulated clusters with varying errors on the temperature prior T e fj. 
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Fig. 1. — l-o" likelihood contours for a simulated region of gas of T e = 10 keV (1.2 x 10 s 
K), v — 200 km/sec, and r = 0.012 obtained by calculating the SZ effect for three different 
observing frequency sets assuming IfiK detector noise. The frequency sets are la) 30, 150, 
300 GHz, lb) 90, 150, 300 GHz, and lc) 145, 225, 265 GHz. The SZ intensity shifts from 
the gas region are calculated using the formula in Nozawa et al. (1998), and the ratios of the 
intensity shifts at different frequencies in each set are computed to eliminate the dependence 
on optical depth. A dot (•) marks the input fiducial gas region. 
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Fig. 2. — Same as Fig. 1 except for v = -200 km/sec. 
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(a) (b) 
Fig. 3. — Same as Fig. 1 except for detector noise of 10/ifT. 
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Fig. 4. — Same as Fig. la except for gas T e and r values: 4a) T e = 7 keV, r = 0.009; 4b) T e 
= 3 keV, r = 0.004. 
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Fig. 5. — Projected 1-cr likelihood contours for the (T e , v, r) parameter space from SZ 
intensity shifts at 145, 225, and 265 GHz for a simulated gas region of T e = 10 keV (1.2 x 10 8 
K), v — 200 km/sec, and r = 0.012 assuming 1/iK detector noise. 
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Fig. 6. — Same as Fig. 5 except for the parameter space (rT e , tv, tT^), corresponding to 
the physical parameters expected to be important. 
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Fig. 7. — Same as Fig. 5 except for the parameter space (a, b, c), the combinations of the 
physical parameters that are best constrained by SZ observations. 





Fig. 8. — SZ simulations of a simulated Nbody+hydro cluster before smoothing and adding 
detector noise. The cluster is about 1O 15 M , has an average gas temperature of about 9 keV, 
and is at z=0.43. Figures 8a, 8b, and 8c are of the 145, 225, and 265 GHz bands respectively. 
Each figure is about 6' x 6' with a pixel size of 0.02' x 0.02'. The images are converted to 
temperature differences from the mean microwave background temperature. The color scale 
is from -lOOfiK to lOOfiK. Primary microwave background fluctuations and point source 
contamination are not included. 
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Fig. 10. — ACT-like SZ simulations of the 9 keV simulated cluster shown in Fig. 8 with 
1' resolution and 3fiK gaussian random instrument noise in each 0.3' x 0.3' pixel. The 
remaining figure specifications are the same as in Fig. 8. 
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(a) (b) (c) 

Fig. 11. — Same as Fig. 10 except using the 3 keV simulated cluster shown in Fig. 9. 




Fig. 12. — The projected 1-a contours of the likelihood surface for the (a e g, b e g, c e s) pa- 
rameter space from the central pixel of simulated ACT-like SZ images of a simulated 9 keV 
Nbody+hydro cluster. The likelihood contours are generated using a Markov chain and are 
for a given noise realization. The star indicates the best fit a c g values from the Markov 
chain. The dot (•) indicates the best fit a e g values obtained with an ideal instrument with- 
out detector noise. A diamond shape (♦) marks the values of line-of-sight integrals through 
the three dimensional cluster. 
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Fig. 13. — The projected l-a contours of the likelihood surface for the (a e ff, b c s, c e s) pa- 
rameter space from the central pixel of simulated ACT- like SZ images of a simulated 3 keV 
Nbody+hydro cluster. The symbols within the contours are as for Fig. 12; note the axes 
are different scales. 




(a) (b) (c) 

Fig. 14. — The projected l-a contours of the likelihood surface for the (a e s, b e s, c e g) param- 
eter space from the central pixel of simulated Planck-like SZ images of a simulated 9 keV 
Nbody+hydro cluster. The symbols within the contours are as for Fig. 12; note the axes 
are different scales. 
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Line-of-Sight Integral 

Fig. 15. — Best fit a e g from simulated, noise-free, ACT-like SZ images found using a Markov 
chain versus J Cdp integrated along the cluster line of sight. The results for four different 
lines of sight through both the 9 keV and 3 keV simulated clusters are plotted. The lines 
of sight are 0', 1', 1.5', and 2' from the central pixel of the simulated SZ images, dp = 
(TdT,vdr,T 2 dr), and C is a matrix of constants introduced in eq. (5). □ = a e g, A = b e g, 
and O = c efr ? and filled (unfilled) shapes correspond to the 9 keV (3 keV) cluster. Typical 
error bars for 1/iK detector noise per 1' beam are ± 0.04, ± 0.3, and ± 3 for a e s, fe e ff, and 
c cS respectively. 



